Reactor z/g disappearance in the Double Chooz experiment 
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The Double Chooz experiment has observed 8,249 candidate electron antineutrino events in 227.93 
live days with 33.71 GW-ton-years (reactor power x detector mass x livetime) exposure using a 
10.3 m^ fiducial volume detector located at 1050 m from the reactor cores of the Chooz nuclear 
power plant in France. The expectation in case of 9i2.= is 8,937 events. The deficit is interpreted 
as evidence of electron antineutrino disappearance. From a rate plus spectral shape analysis we find 
sin^ 2^13 = 0.109 ± 0.030(stat) ± 0.025(syst). The data exclude the no-oscillation hypothesis at 
99.8% CL (2.9cr). 



I. INTRODUCTION 

In the three neutrino paradigm, there are three mix- 
ing angles that can be measured in neutrino oscillation 
experiments. For many years, the CHOOZ reactor neu- 
trino experiment [l| had the best hmit on the value of 
013. Recently, the value of ^13 has been shown to be 
non-zero by the combination of fits to KamLAND and 
solar [13, MINOS Q, T2K @ and, more precisely, by 
the new generation of reactor antineutrino disappearance 
experiments: Double Chooz [7], Daya Bay [8] and RENO 

The Double Chooz analysis is unique among reactor 
experiments in its fit to the energy spectrum. In the pre- 
vious reactor measurements of ^13, Double Chooz pre- 
sented both a rate-only analysis and an analysis using 
both the rate and the shape of the energy spectrum, 
while Daya Bay and RENO presented rate-only analyses. 
The disappearance of reactor electron antineutrinos has 
a well-defined effect on the shape of that spectrum. The 
use of the energy distribution to constrain the oscillation 
parameters requires a good understanding of the energy 
response of the detector and of the accuracy of the Monte 
Carlo. That understanding is achieved through multiple 
calibration techniques, in time, space and energy. 

This paper continues the analysis reported in Q with 
a larger data set, a new energy scale definition, reduced 
background rates and improved systematic uncertainties. 
Additionally, the running period has been subdivided 
into a two-reactor-on period and a one-reactor-on period 
in the oscillation fit to help separate signal and back- 
ground. 

Reactor antineutrinos are observed using the inverse 
beta decay (IBD) reaction z/e + p — > e+ + n in which 
there is a positron whose signal is promptly seen, and a 
neutron, whose delayed signal is seen after a mean time of 
about 30 /is from its capture in the gadolinium-doped tar- 
get. The prompt energy of the positron allows us to de- 
termine the antineutrino energy and observe the antineu- 
trino spectrum. The energy deposited by the positron 
including annihilation is related to antineutrino energy 
E^^ by -Eprompt — Ejj^ — T„ — 0.8 MeV where Tn denotes 
the average neutron recoil energy and is small compared 

to Ey^ . 

The previous analysis represented 15.34 GW-ton-years 
of exposure, taking into account the reactor livetime and 



the detector fiducial mass. Here we re-analyze that data 
set together with an additional 18.37 GW-ton-years giv- 
ing a total of 33.71 GW-ton-years. In addition the anal- 
ysis of 22.5 hours of both-reactors-off data allows a cross 
check of our estimates of the correlated and accidental 
backgrounds. 

The structure of the paper is as follows. In Section 
iniwe review the experimental setup and detector. Sec- 
tion mil covers the measurements and simulations of the 
Chooz reactors used to predict the unoscillated neutrino 
spectrum, as well as the model used to describe the de- 
tector. Event reconstruction including the energy deter- 
mination of candidate events is described in Section llVI 
The steps that are used to identify reactor neutrino can- 
didates are covered in Section [V] Section IVTl presents the 
extraction of neutrino mixing parameters from the mea- 
sured antineutrino rate and energy distribution. 

II. DETECTOR AND METHOD DESCRIPTION 
A. Overview 

The Double Chooz detector system [l^ consists of a 
main detector, an outer veto, and calibration devices 
(Figure [T]). The main detector is made of four concentric 
cylindrical tanks with a chimney in the center at the top 
and is filled with liquid scintillators or mineral oil. 

The innermost 8 mm thick transparent (UV to visible) 
acrylic vessel contains 10.3 m'^ gadolinium loaded liquid 
scintillator called the (/-target (NT). The NT volume is 
surrounded by the 7-catcher (GC), a 55 cm thick Gd-free 
liquid scintillator layer in a second 12 mm thick acrylic 
vessel, used to detect gamma rays escaping from the v- 
target. Outside the 7-catcher is the buffer, a 105 cm thick 
mineral oil layer. It shields from radioactivity of photo- 
multiplier tubes (PMTs) and surrounding rock, and is 
one of the major improvements over the CHOOZ exper- 
iment [l|. The 390 10-inch PMTs [IMll are installed 
on the inner wall of the stainless steel buffer tank to col- 
lect light from the inner volumes. These three volumes 
and PMTs constitute central detector system referred to 
as the inner detector (ID). Outside the ID, and optically 
separated from it by a stainless steel vessel, is a 50 cm 
thick inner veto (IV) liquid scintillator. It is equipped 
with 78 8-inch PMTs and functions as a cosmic muon 
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veto and as an active shield to spallation neutrons pro- 
duced outside the detector. The detector is covered and 
surrounded by 15 cm of demagnetized steel to suppress 
external gamma rays. The main detector is covered by 
an outer veto system (OV) described in Section Hi Gl 



B. Radiopurity 

All parts of the Double Chooz detector have been 
thoroughly screened for their content of radioactive iso- 
topes prior to their installation. The screening was 
carried out by direct gamma spectroscopy with a va- 
riety of germanium detectors in underground laborato- 
ries. Among them were the large HPGe detector for the 
non-destructive radioassay at Saclay and the GeMPI 
detector at Gran Sasso [la] with a sensitivity of about 
10 fiBq/kg for U and Th. In addition, neutron activation 
analyses have been performed for dedicated parts of the 
inner detector: the acrylics for NT and GC vessels as well 
as the wavelength shifter PPG 16] . The irradiations were 
done at the FRM II research reactor in Garching, Ger- 
many by a thermal neutron flux of 1.63T0^^ cm~^-s~^, 
with subsequent gamma spectroscopy in the Garching 
underground laboratory [l7| . 

The PMT glass and cavern rock are the main sources 
of the gamma ray background. The PMT glass was made 
from low activity sands using a platinum coated furnace 
to reduce contamination. Radioactivity of the glass sam- 
ples was measured during development of the low activ- 
ity glass and production of the PMTs [l^ . The average 
measurements were 13 ppb, 61 ppb and 3.3 ppb for ^^®U, 
^■^^Th and *'^K, respectively assuming radio-equilibrium, 
which are much smaller than regular PMT glass. 
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FIG. 1. A cross-sectional view of the Double Chooz detector 
system. 




The design goal of Double Chooz concerning radiopu- 
rity is no more than ^0.8 accidental background events 
per day. Along with the radiopurity screenings. Dou- 
ble Chooz maintained strict clean-room conditions dur- 
ing the setup of the detector with an ISO-level up to 6. 
The analysis of BiPo coincidences in the detector data 
yields concentrations of U and Th in NT and GC be- 
low the design goal of 10~^^g/g. The accidentals rate is 
measured to be < 0.5 d~^, well below our design goal. 
The daily rate of correlated background events stemming 
from (a, n)-reactions of ^^°Po on ^^C is estimated to be 
smaller than 0.020 d~^ (scaled from the result of Kam- 
LAND [11]), which is negligibly small compared to the 
neutrino signal. 



C. Double Chooz Liquids 

The CHOOZ experiment was limited in sensitivity by 
the optical instability of its gadolinium-loaded (Gd) scin- 
tillator Therefore a new type of metal loaded or- 
ganic liquid scintillator was developed for Double Chooz 
[la |. The target scintillator used in the NT must ful- 
fill the basic requirements of Gd solubility in the solvent 
of choice, optical transparency, radiopurity and chemical 
stability. In addition, the organic liquid must be com- 
patible with the detector materials in contact with the 
scintillator, mainly acrylics. Safety considerations influ- 
enced the scintillator design as well. 

Since the rare earth Gd does not dissolve in the re- 
quired amount in the organic solvents used for liquid 
scintillators, a metalorganic complex is formed provid- 
ing higher solubility. In particular, the complex of choice 
is a metal- ;3-diketone, Gd(thd)3, Gd(III)-tris-(2,2,6,6- 
tetramethyl-heptane-3,5-dionate). Such complexes are 
known for their stability and high vapor pressure. This 
allowed us to purify the material by sublimation reducing 
radioimpurities U, Th and K. The Gd concentration in 
the NT is 0.123% by weight, which corresponds to about 
1 g/liter. 

As scintillator solvent for the NT we have chosen 
an ortho-phenylxylylethane (o-PXE)/n-dodecane mix- 
ture at a volume ratio of 20/80. To shift the scintil- 
lation light into a more transparent region, wavelength 
shifters are added. In both scintillators we use PPO (2,5- 
diphcnyloxazole) as primary fluor and bis-MSB (4-bis-(2- 
mcthylstyryl)benzene) as secondary wavelength shifter. 

The hght yield and density of the GC hquid (22.5 m^) 
were matched simultaneously to the NT values ^2)]. To 
achieve this goal, a medicinal white oil was added as a 
third solvent to the GC. The light yield of the GC is op- 
timized for homogeneous detector response using Monte 
Carlo simulations. To avoid mechanical stress on the de- 
tector vessels the densities of all four detector liquids were 
matched at the detector temperature of about 15 °C to 
0.804 ± 0.001 g/cm^. 

The attenuation lengths for wavelengths in the region 
of scintillator emission are well above the dimensions of 
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FIG. 2. Average target detector response evolution in time, 
as measured by the mean energy of the Gd-capture peak aris- 
ing from interaction of spaUation neutrons in the NT. 



the corresponding vessels. Optical stability of the scin- 
tillators is demonstrated in Figure [2l where the stability 
of the peak energy of neutron captures on Gd is shovifn. 
The energy response of the detector was found to be sta- 
ble within 1% over the data-taking period of about one 
year. 

The absolute number of H nuclei ("proton number") 
as well as the precision on its knowledge are crucial pa- 
rameters. The error on the proton number is minimized 
by using well defined and pure chemicals in combination 
with a precise knowledge of the weights of each chemi- 
cal added in the scintillator production. The amount of 
NT scintillator was determined after thermalization by 
a weight measurement with a precision of 0.04%. The 
hydrogen fraction in the NT is 13.6% by weight, known 
with 0.3% relative precision. This error includes the un- 
certainties originating from the weights of the scintillator 
ingredients. In addition, the error takes into account the 
knowledge of the hydrogen content of not fully defined 
impurities in the chemicals which are on the per mil level 
for the main components [T6|. 

A mixture of solvents was used in all detector volumes 
to allow for density matching. The 110 m'^ of buffer liq- 
uid contain a medicinal white oil (53% by volume) and 
an n-alkane mixture (47% by volume). This liquid was 
optimized for transparency and low aromaticity to mini- 
mize scintillation light production in the buffer. The veto 
volume is filled with 90 m'^ of liquid scintillator, a mix- 
ture of linear alkyl benzene (LAB) and n-alkanes, with 
2 g/1 PPO as fluor and 20 mg/1 bis-MSB as secondary 
wavelength shifter. 



D. ID Photomultiplier Tubes 

The inner detector uses 390 Hamamatsu R7081 10- inch 
PMTs [13 to view the target volume. The glass is a low 
background type, contributing only a few Hz of singles 
rate in the detector. The PMTs are operated with a 



gam of 10'' at the PMT anode. They are submerged in 
a paraffin oil buffer liquid. The base circuit is enclosed 
in a transparent epoxy resin. Some PMTs are observed 
to generate light flashes from their base circuit through 
the epoxy resin, causing false triggers. HV for the 14 
worst PMTs was turned off. Since the signal pattern 
is different from that of the neutrino signal, the false 
events are safely removed from the neutrino sample as 
described in Section HVP] The 800 PMTs for both this 
and an eventual near detector were characterized care- 
fully [Hill, [H. The following characteristics were mea- 
sured: for one photoelectron signals, the ratio of the one 
photoelectron peak to the valley between that peak and 
the pedestal was 4, with 1/4 photoelectron thresholds; 
the quantum efficiency x collection efficiency (efficiency 
that photoelectrons produced in the cathode are collected 
by the first dynode) was 23%; transit time spread was 
3 ns (FWHM); the afterpulse probability was in average 
2.7%; the charge output was linear up to 300 photoelec- 
trons per PMT; dark hit rate was approximately 2 kHz 
measured 20 hours after turning on the HV. Each PMT is 
shielded by a mu-metal cylinder to suppress effects from 
the gamma shield and the earth's magnetic field [l3| and 
is equipped with an angle- adjust able mounting jig. The 
PMTs are angled to collect light more uniformly from the 
detector. 



E. The Inner Veto 

The IV is a cylindrical stainless steel vessel (radius 
3.3 m and height 6.8 m) surrounding the ID and optically 
separated by the buffer tank. It shields the ID with a 
50 cm thick layer of liquid scintillator against external 
radioactivity and spallation neutrons created by cosmic 
muons. At the same time it acts as an active detec- 
tor identifying cosmic muons crossing it. The design of 
the IV was optimized by the use of a MC simulation 
24| . where the emphasis was on a high number of de- 
tected photoelectrons (PE) per MeV deposited in the IV 
volume and on a high efficiency in rejecting muons and 
correlated background events produced by them. The 
resulting configuration of the IV consists of 78 PMTs, 
divided into three parts: the top has 24 PMTs, the side 
walls have 12 PMTs at the mid way point and the bot- 
tom has 42 PMTs. The 78 8-inch PMTs (Hamamatsu 
R 1408), which were previously used in the 1MB and 
Super-Kamiokande experiments, were tested and modi- 
fied for use in Double Chooz [2?|. Each IV PMT and 
its base are contained in a stainless steel encapsulation, 
with a transparent PET window at the front end. The 
capsules are filled with mineral oil to match the opti- 
cal properties of the surrounding scintillator. All sur- 
faces of the IV are painted with highly reflective white 
coating (ARIOO/CLX coating from MaxPerles [26j), the 
side walls of the buffer vessel are covered with reflective 
VM2000 sheets. Using the OV, the muon rejection ef- 
ficiency was found to be larger than 99.99% for muons 
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crossing the IV volume. 



F. Electronics and Data Acquisition 

The full readout and data acquisition (DAQ) for both 
the ID and the IV detectors are depicted in Figure [31 
The functional principle is that digitization of PMT sig- 
nals (see Section III Pp is done by flash- ADC electronics. 
As shown in Figure [3l from left to right, the electron- 
ics elements are the High Voltage (HV) splitter, the HV 
supply, the Front- End electronics (FEE), the Trigger sys- 
tem [27] and the flash- ADC digitizing electronics [2^, [20| 
(i^-FADC). Each PMT has a single cable for both PMT 
signal (5 mV per PE) and HV (~ 1.3 kV). A custom made 
HV-splitter circuit decouples both components. The HV 
is provided by CAEN-A1535P [30] supplies. PMT sig- 
nals are optimized (amplified, clipped, baseline restored 
and coherent noise filtered) by the FEE for digitization. 
The FEE also delivers sum signals, whose amplitude is 
proportional to charge, that are fed into a custom trig- 
ger system. The circuit generating the sum signal sub- 
tracts the input amplitude after about 100 ns. This ca- 
pability allows the trigger input signals to suffer from 
minimal overshoot that can lead to trigger dead-time. 
This same feature works as a high-pass filter: slow sig- 
nals (frequency ^ 1 MHz) cannot cause a trigger. The 
ID PMTs are separated into two ID super-groups at the 
trigger level, uniformly distributed across the volume. 
Either super-group can cause a trigger of the ID based 
on energy and sub-group multiplicity information. The 
ID triggers at energies about 350 keV. The trigger effi- 
ciency is 100.0% above the analysis threshold 0.7 MeV 
with negligible uncertainty. Both energy and sub-group 
multiplicity information are used to cause IV triggers. 
The IV triggers at ^ 10 MeV which corresponds to 8 cm 
of a minimum ionizing muon track. The v-FADG sys- 
tem relies on 64 CAEN-Vxl721(VME64x) fsd] waveform 
digitizers. Each card has 8 channel with 8-bit flash- ADC 
(FADC) at 500 MS/s. Each channel holds up to 1024 
4 fj,s waveforms without readout. When triggered, the 
256 ns waveform is recorded, containing > 90% of the 
scintillation light emitted. Up to ~ 3 MeV, a single-PE 
is deposited per channel, each having ~ 40 mV amplitude 
corresponding to around 10 samples per PE. FADC am- 
plitude saturation leads to some degree of non-linearity 
for > 100 MeV energies. Above 500 MeV, up to a 40% 
non-linearity has been estimated. 

The FADC baselines are observed to be stable, show- 
ing variations below 1 ADC. After power-cycling, small 
(sub-mV) DC shifts in the baseline are observed. Due 
to under-sampling of the baseline, these shifts can cause 
a bias in the reconstructed charge estimation. This bias 
manifests itself as an effective non-linearity for signals 
below 2 PEs and has been thoroughly studied, measured 
and calibrated out, as described in Section Section HV El 

All systems (trigger, ID and IV) are readout by the 
same DAQ upon any trigger of either the ID or IV. The 
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FIG. 3. Block diagram of the Double Chooz readout and 
DAQ systems. 



system is deadtime free, as demonstrated by two monitor 
systems running at 2 Hz. The dead time monitor wave- 
forms are, in addition, used to randomly sample the de- 
tector providing extra baseline monitoring, background 
and dark-current information. 



G. The Outer Veto 



The OV is installed above the ID, IV and 15 cm of 
shielding steel. A lower outer veto is mounted directly 
above the shielding and provides (z, y) coordinate for 
muons passing through a 13 m x 7 m area centered 
on the chimney; a 110 cm x 30 cm region around the 
chimney is left open. The lower outer veto has been in- 
stalled for 68.9% of the data presented here, and is used 
to help reduce background levels quoted in Q. An up- 
per outer veto, again measuring (a;, y) coordinates, has 
been mounted above the chimney and glove box used for 
source insertion, to cover this area. The upper outer veto 
was not present for this analysis. 

The outer veto is assembled from modules containing 
64 scintillator strips, each 5 cm x 1 cm x 320 cm or 360 
cm. Each strip was extruded with a hole running through 
its length, through which a 1.5 mm diameter wavelength- 
shifting fiber was threaded. Modules are built out of two 
superimposed 32-strip layers with the top layer offset by 
2.5 cm from the bottom layer. The 64 fibers are coupled 
at one end to a Hamamatsu H8804 multi-anode photo- 
multiplier tube (M64); the other fiber ends are mirrored. 
The OV modules are positioned over the inner detector 
in two layers, one with strips oriented in the x direction 
and one in the y direction. Each M64 is connected to a 
custom front-end board with a MAR0C2 ASIC [sH and 
an FPGA. The MAR0C2 allows adjustment of the elec- 
tronic gain of each of the 64 channels, which is needed 
to correct for the factor of 2 pixel-to-pixel gain variation 
in the M64. Signals that exceed a common threshold 
are sent to a multiplexed 12-bit ADC, providing charge 
information for hit strips. 
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H. Calibration Systems 

The experiment is calibrated using light sources, ra- 
dioactive (point-like) sources, and cosmic rays. 

A multi-wavelength LED-fiber system (LI) is used to 
inject light into the inner detector and the inner veto 
from a set of fixed points. The optical fibers of the LI 
are routed inside the detector and the fiber ends are at- 
tached to the PMT covers. Some of the injection points 
are equipped with diffuser plates to widely illuminate the 
detector. The other ends of the fibers are connected to 
blue and UV LEDs (385, 425 and 470 nm wavelengths 
for the ID, and 365 and 475 nm for the IV) whose flash 
rate, light intensity and pulse width are controlled re- 
motely. Data are taken with the LI systems regularly. 
The LI data are used to measure the PMT and readout 
electronics gains and the time offsets and to monitor the 
stability of those gains and offsets. 

Radio-isotopes ^^"^Cs, ^^Qe, ^oco, and '^^'^Ci, sealed 
in miniature capsules, have been deployed in the NT 
and GC. The visible energy response is measured with a 
0.662 MeV gamma (Cs-137), 2x0.511 MeV annihilation 
gammas (Ge-68), which also corresponds to the thresh- 
old for inverse beta decay, the combination of 1.173 MeV 
and 1.333 MeV gammas (Co-60), and the 2.223 MeV 
gamma from neutron capture on hydrogen (Cf-252). The 
detector response to neutrons is calibrated using ^^^Cf. 
Source rates are at the level of 50 Bq. 

Deployments in the NT are realized by lowering the 
sources from a glove box at the detector top through 
the detector chimney. A motorized puUey-and-weight 
system, operated from a glove box, is used to position 
sources at positions along the target symmetry axis. The 
range of deployments is from 1 cm above the NT bottom 
up to the chimney; the positions of the source are known 
within 1 mm. In the GC, the source is attached to a 
motor-driven wire and guided through a rigid hermetic 
looped tube (GT). The sources are inserted in the GT 
near the chimney top. The loop traverses interior re- 
gions of the GC and passes near boundaries with the NT 
and the buffer. The position of the source along the loop 
is known to 1 cm, and in the NT boundary region, the 
perpendicular distance between the source and the target 
wall is known within 2 mm. The materials of the source 
capsules and deployment systems in the NT and GC are 
modeled by the detector simulation. 

Cosmic rays are analyzed to identify stopping muons, 
spallation neutrons, and cosmogenic radioactive isotopes. 
Several thousand spallation neutrons per day are cap- 
tured on hydrogen and gadolinium in the ID. 

The use of the calibration data for issues of energy 
uniformity, stability, non-linearity and absolute calibra- 
tion is described in Section flV El The neutron detection 
efficiency from ^^^Cf is described in Section FV Fl Good 
control of uncertainties on detection efficiency is essen- 
tial for sensitivity to neutrino disappearance with a single 
detector. The detailed calibration data allow a precise 
energy-shape fit to the prompt neutrino candidates for 



the most sensitive extraction of ^13. 



III. REACTOR AND DETECTOR MODELS 



A. Thermal Power 



Double Chooz's sources of antineutrinos are the reactor 
cores Bl and B2 at the Electricite de France (EDF) Cen- 
trale Nucleaire de Chooz. Antineutrinos are produced 
in nuclear reactors by the /3-decay of the fission prod- 
ucts. Four main isotopes, ^^^U, ^^^Pu, ^^*U, and ^^^Pu, 
provide >99.7% of the fissions and antineutrinos. 

Chooz Bl and B2 are N4 type pressurized water re- 
actor (PWR) cores, and as such are two of the most 
powerful cores in the world with nominal thermal power 
outputs of 4.25 GWth each. The instantaneous thermal 
power of each reactor core is provided by EDF as 
a fraction of the total power and is evaluated over time 
steps of <1 minute. The instantaneous thermal power is 
derived from the in-core instrumentation with the most 
important variable being the temperature of the water in 
the primary loop. 

The in-core instrumentation calibration is tested 
weekly using the heat balance in the secondary loop, 
which is heated by the primary loop containing water 
heated by fissions. In the secondary loop, steam is gen- 
erated to drive turbines. By using measurements of the 
heat flow in the secondary loop, the thermal power can 
be measured. This test is performed with the reactor 
running at full power. The uncertainty at lower power is 
therefore slightly larger. The in-core instrumentation is 
re-calibrated if it deviates by more than the uncertainty 
in the heat balance measurement. 

Since the accuracy of the thermal power measurement 
determines the maximum power at which the core can 
operate, EDF has performed a detailed study of the un- 
certainty in this measurement (32|-|3^. The dominant 
uncertainty on the weekly heat balance at the secondary 
loops comes from the measurement of the water flow. At 
the nominal full power of 4250 MW the final uncertainty 
is 0.5% (1 cr C.L.). Since the amount of data taken with 
one or two cores at intermediate power is small, this un- 
certainty is used for the mean power of both cores. This 
is smaller than the typical uncertainty for PWRs of 0.7% 
[3^ and reflects optimizations in the pipe geometry of the 
secondary loop, as well as great care taken to understand 
the sensor uncertainties, including full-scale test stands 
for the most critical sensors. 
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B. Mean Cross Section per Fission 

The mean cross section per fission is effectively a spec- 
trum averaged cross section. It is given by 

/"OO 

(^/) = y]"'^^'^/'''^ ^ / dE Sk{E)(JiBD{E) 

k k ^0 

(1) 

where is the fractional fission rate of the fc"* isotope 
(k = 235u^ 239py^ 238u^ 241 p^)^ ^j^g reference 

spectrum of the fc*'' isotope and amo is the inverse beta 
decay cross section. The determinations of the ak require 
the simulation of the reactor core fSection llll C|) . 

The antineutrino spectrum for each fission isotope is 
the result of the beta decays of many different fission 
products. For ^^^U, ^^^Pu, and ^''^Pu, the reference an- 
tineutrino spectra are derived from measurements of the 
/3 spectra at the ILL research reactor [36l-[38j. In the 
case of an ab initio calculation of the spectrum is 

used |3§|. The conversion of the /3 spectra to antineu- 
trino spectra has recently been improved by using more 
data on the many l3 transitions and higher order energy 
corrections [H, \^ . We use the conversion scheme of (ioj 
including corrections for off-equilibrium effects [4l|. The 
uncertainty on these spectra is energy dependent but is 
on the order of 3%. The new technique for the analysis of 
the /3 spectra has led to an overall change in the normal- 
ization of the Sk{E) that, when applied to previous re- 
actor antineutrino experiments, results in measurements 
that are lower than predictions for experiments at short 
baselines [41,] . 

C. Fission Rate Computation 

The fractional fission rates ak of each isotope are 
needed in order to calculate the mean cross section per 
fission of (Equation [T|). They are also required for the 
calculation of the mean energy released per fission for 
reactor R: 

{Ef)R^Y.'^k{Ef)k. (2) 

k 

The mean energies released per fission per isotope {Ef)k 
are summarized in Table HI The thermal power one would 
calculate given a fission is relatively insensitive to the 
specific fuel composition since the {Ef)k differ by <6%; 
however, the difference in the detected number of an- 
tineutrinos is amplified by the dependence of the norm 
and mean energy of S'fe(E) on the fissioning isotope. For 
this reason, much effort has been expended in developing 
simulations of the reactor cores to accurately model the 
evolution of the . 

Double Chooz has chosen two complementary codes 
for modeling of the reactor cores: MURE and DRAGON 
[13-01 • MURE is a 3D full core simulation which uses 
Monte Carlo techniques to model the neutron transport 



in the core. DRAGON is a 2D simulation which models 
the individual fuel assemblies. Using some approxima- 
tions, it solves the neutron transport equation in the core. 
These two codes provide the needed flexibility to extract 
fission rates and their uncertainties. These codes were 
benchmarked against data from the Takahama-3 reactor 
and were found to be consistent other codes commonly 
used in the reactor industry for reactor modeling within 
the uncertainty in the Takahama data [46[. 

The construction of the reactor model requires detailed 
information on the geometry and materials comprising 
the core. The Chooz cores are comprised of 205 fuel 
assemblies. For every reactor fuel cycle, approximately 
one year in duration, one third of the assemblies are re- 
placed with assemblies containing fresh fuel. The other 
two thirds of the assemblies are redistributed to obtain a 
homogeneous neutron flux across the core. The Chooz re- 
actor cores contain four assembly types that differ mainly 
in their initial 235y enrichment. These enrichments are 
1.8%, 3.4% and 4%. 

The data set presented here spans fuel cycle 12 for core 
B2 and cycle 12 and the beginning of cycle 13 for Bl. 
EDF provides Double Chooz with the locations and ini- 
tial burnup of each assembly. Based on these maps, a full 
core simulation was constructed using MURE for each cy- 
cle. In addition, the beginning-of-fuel-cycle composition 
needs to be determined based on the burnup of each as- 
sembly. To accomplish this, an assembly-level reference 
simulation is run using both MURE and DRAGON for 
each of the four fuel assembly types. The results of the 
reference simulations are compared to EDF's own simula- 
tion code AP0LL02-F from which the burnup values are 
derived. The uncertainty due to the simulation technique 
is evaluated by comparing the DRAGON and MURE re- 
sults for the reference simulation leading to a small 0.2% 
systematic uncertainty in the fission rate fractions a^. 

Once the initial fuel composition of the assemblies is 
known, MURE is used to model the evolution of the full 
core in time steps of 6 to 48 hours, depending on the op- 
erating conditions of the reactor. The results from each 
simulation time step are written to a database. This al- 
lows the afc's, and therefore the predicted antineutrino 
flux, to be calculated. The results averaged over the cur- 
rent data set are shown in Table [J 

The systematic uncertainties on the a^'s are deter- 
mined by varying the inputs and observing their effect 
on the fission rate relative to the nominal simulation. 
The uncertainties considered are those due to the thermal 
power, boron concentration, moderator temperature and 
density, initial burnup error, control rod positions, choice 
of nuclear databases, choice of the energies released per 
fission, and statistical error of the MURE Monte Carlo. 
The systematic errors associated with each input are con- 
sidered independently and the uncertainties propagated 
quadratically. The correlation coefficients among isotopic 
fission rates due to the thermal power constraint are also 
computed, and a covariance matrix is constructed with 
these contributions in order to properly account for those 
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TABLE I. Mean energy released per fission {Ef)k from [431 
and fractional fission rate {ak) of the isotope k for this data. 



Isotope {Ef)k (MeV) 

— 2351 



201.92±0.46 
209.99±0.60 
205.52±0.96 
213.60±0.65 



0.496±0.016 
0.351±0.013 
0.087±0.006 
0.066±0.007 



correlations. The uncertainties in the a^'s are listed in 
Table HI The two largest contributions come from the 
moderator density and control rod positions. 



TABLE IL The uncertainties in the antineutrino prediction. 
All uncertainties are assumed to be correlated between the 
two reactor cores. They are assumed to be normalization and 
energy (rate and shape) unless noted as normalization only. 



Source 


Normalization Only 


Uncertainty [%] 


Pth 


yes 


0.5 




yes 


1.4 




no 


0.2 


{Ef) 


no 


0.2 


Lr 


yes 


<0.1 


„ R 


no 


0.9 


Total 




1.8 



D. Bugey4 Normalization and Antineutrino Rate 
Calculation 



In the current, far-only, phase of Double Chooz, the 
rather large uncertainties in the reference spectra of Sec- 
tion |TirB] limited our sensitivity to Bi^. To mitigate this 
effect, the normalization of the cross section per fission 
for each reactor is "anchored" to the Bugey4 rate mea- 
surement at 15 m 



Bugey 



){<yf)k 



(3) 



where R stands for each reactor. The second term 
corrects for the difference in fuel composition between 
Bugey4 and each of the Chooz cores. This treatment 
takes advantage of the high accuracy of the Bugey4 an- 
chor point (1.4%) and suppresses the dependence on the 
predicted {(Jf)R- This is due to the smallness of the cor- 
rection term (a^ — a^"^*^^). At the same time, the analy- 
sis becomes insensitive to possible oscillations at shorter 
baselines due to heavy Am^ ~ 1 eV^ sterile neutrinos. 

The expected number of antineutrinos with no oscilla- 
tion in the i^^ energy bin with the Bugey4 anchor point 
becomes: 



eNr. 1 



pR 



4^ {Ef)i 



(4) 



where e is the detection efficiency, Np is the number of 
protons in the target, Ljj is the distance to the center of 
each reactor, and is the thermal power. The variable 
{Ef)R is the mean energy released per fission defined in 
Equation [21 while {(Jf)^ is the mean cross section per 
fission defined in Equation [3l The three variables P^, 
{Ef)R and {<Jf)R are time dependent with {Ef)R, and 
{(Tf)R depending on the evolution of the fuel composition 
in the reactor and Pff^ depending on the operation of the 
reactor. 

A covariance matrix Mff^ = SNf^^SNf^^ is con- 
structed using the uncertainties listed in Table [III This 



matrix is constructed in terms of real energy and is 
converted into reconstructed energy by running multi- 
ple simulations drawn from a Cholesky decomposition of 



exp 



For these simulations, the full detector Monte 



Carlo described below is used. The use of Equation [4] to 
construct the covariance matrix allows time and spectral 
information to propagate to the final analysis. 

The IBD cross section used is the simplified form from 
Vogel and Beacom 49] : 



mt 



(5) 



where 



Ee+ — — 




Anip -E^ + A - 



A2 



2m„ 



(6) 

and rrig and i?e+ are the positron mass and energy. The 
variables m„ and nip are the masses of the neutron and 
proton with A = to„ — rup. The constant K is in- 
versely proportional to the neutron lifetime. We use the 
MAMBO-II measurement of the neutron lifetime [10] and 
find K = 0.961 X 10-^3 ,,^^2 MeV'^. 



E. Detector Model 

We model the detector response using a detailed 
Geant4 [sij simulation with enhancements to the scintil- 
lation process, photocathode optical surface model, and 
thermal neutron model. Apart from these additions, the 
physics list is similar to Geant4's QGSP_BERT_HP refer- 
ence physics list [s^ . without processes for high- mass 
hadrons. Our custom scintillation process implements 
detailed light waveforms, spectra, re-emission, and Birks- 
law [s^ quenching. Our photocathode model is based on 
a standard mathematical model of a thin, semitranspar- 
ent surface with absorption and refractive index (s^ , and 
also includes the collection efficiency for photoelectrons 
as a function of position of emission on the photocathode. 
Our custom neutron thermalization process implements 
molecular elastic scattering for neutrons under 4 eV and a 
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radiative capture model with improved final state gamma 
modeling. 

The simulation models the detector geometry to a fine 
level of detail, particularly with regard to the geome- 
try of the phototubes and mu-metal shields and of all 
materials near the active volume such as tank walls and 
supports. The orientation and positions of the phototube 
assemblies were set using data from a photographic sur- 
vey with sub-mm accuracy. The dimensions of the tank 
walls and supports were checked by experimenters during 
assembly and installation, and placement also verified by 
photographic survey. 

Simulated IBD events are generated with run-by-run 
correspondence of MC to data, with fluxes and rates cal- 
culated as described in Section UlIDI Radioactive decays 
in calibration sources and spallation products were simu- 
lated using detailed models of nuclear levels, taking into 
account branching ratios and correct spectra for transi- 
tions [55]. 

Optical parameters used in the detector model are 
based on detailed measurements made by the collabo- 
ration. The relative light yield of the NT compared to 
the GC was measured using a Compton backscatter peak 
method in order to select scattered electrons with fixed 
energy [2l|. Tuning of the absolute and relative light 



yield in the simulation was done with calibration data. 
The scintillator emission spectrum was measured using 
a Gary Eclipse fluorometer [s^. The photon emission 
time probabilities used in the simulation are obtained 
with a dedicated laboratory setup [16]. For the ioniza- 
tion quenching treatment in our MG, the light output of 
the scintillators after excitation by electrons [57] and al- 
pha particles [58i] of different energies was measured. The 
non-linearity in light production in the simulation has 
been adjusted to match these data. The attenuation and 
re-emission probabilities of each of the scintillator compo- 
nents in the relevant wavelength range are implemented 
in the MG. The fine-tuning of the total attenuation was 
made using measurements of the complete scintillators 
[l6| . Other measured optical properties include reflectiv- 
ities of various detector surfaces and indices of refraction 
of detector materials. 



F. Readout System Simulation 

The Readout System Simulation (RoSS) accounts for 
the response of elements associated with detector read- 
out, such as from the PMTs, FEE, FADGs, trigger sys- 
tem and DAQ. The simulation relies on the measured 
probability distribution function (PDF) to empirically 
characterize the response to each single PE as measured 
by the full readout channel. The Geant4-based simu- 
lation calculates the time at which each PE strikes the 
photocathode of each PMT. RoSS converts this time-per- 
PE into an equivalent waveform as digitized by FADGs. 
A dedicated setup was built to measure most of the nec- 
essary PDFs as well as to tune the design of the full 



readout chain. Ghannel-to-channel variations, such as 
gains, baselines, noise, single PE widths, etc., are taken 
into consideration, to accurately predict dispersion ef- 
fects. This capability allows the simulation to exhibit 
non-linearity effects as observed in the data, as described 
in Section IIIFI After calibration, the MG and data ener- 
gies agree within 1%. About 25% of the width of the cal- 
ibrated H-capture (2.2 MeV 7 line) results from readout 
effects; i.e., effects beyond photon-statistics fluctuations. 



G. Monte Carlo Event Generation 

A set of Monte Garlo i>e events representing the ex- 
pected signal for the duration of physics data-taking is 
created based on the formalism of Equation 01 The cal- 
culated IBD rate is used to determine the rate of inter- 
actions. Parent fuel nuclide and neutrino energies are 
sampled from the calculated neutrino production ratios 
and corresponding spectra, yielding a properly normal- 
ized set of IBD-progenitor neutrinos. 

Once generated, each event-progenitor neutrino is as- 
signed a random creation point within the originating re- 
actor core. The event is assigned a weighted-random in- 
teraction point within the detector based on proton den- 
sity maps of the detector materials. In the center-of-mass 
frame of the v — p interaction, a random positron direc- 
tion is chosen, with the positron and neutron of the IBD 
event given appropriate momenta based on the neutrino 
energy and decay kinematics. These kinematic values 
are then boosted into the laboratory frame. The result- 
ing positron and neutron momenta and originating ver- 
tex are then available as inputs to the Geant4 detector 
simulation. "Truth" information regarding the neutrino 
origin, baseline, and energy are propagated along with 
the event, for use later in the oscillation analysis. 



IV. RECONSTRUCTION 
A. Pulse Reconstruction 

The pulse reconstruction provides the signal charge 
and time in each PMT. Pulser triggers are taken with 
a rate of 1 Hz in order to provide accurate information 
about the baseline for each of the 468 readout channels. 
The baseline mean (i?mean) and rms (-Brms) are com- 
puted using the full readout window (256 ns). 

The integrated charge (g) is defined as the sum of dig- 
ital counts in each waveform sample over the integration 
window, once the pedestal has been subtracted. The 
pedestal is computed as the integration of -Bmcan over 
the same window. In order to improve the charge reso- 
lution, the size of the integral window has been set to a 
112 ns subsample of the readout one, based on the width 
of the single PE signals. 

In order to find the pulses within the readout win- 
dow, a dynamic window algorithm is used. The algo- 
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rithm searches for the 112 ns window which maximizes 
the intcgraL In the absence of an actual PE signal, this 
algorithm would reconstruct the largest noise fluctuation, 
leading to a bias in the charge reconstruction. To address 
this, we introduce two requirements: > 2 ADC counts in 
the maximum bin, and q > Brms x v'^; where Ng is the 
number of integrated waveform samples (56 for a 112 ns 
window). For each pulse reconstructed, the start time is 
computed as the time when the pulse reaches 20% of its 
maximum. This time is then corrected by the PMT-to- 
PMT offsets obtained with the LI system. 



B. Vertex Reconstruction 

Vertex reconstruction in Double Chooz is not used for 
event selection, but is used for event energy reconstruc- 
tion. It is based on a maximum charge and time likeli- 
hood algorithm which utilizes all hit and no-hit informa- 
tion in the detector. Assuming the event to be a point- 
like source of light characterized by the set 

^ = {xo,yo,zo,to,^) (7) 

where (xq, yo,zo) is the event position in the detector, to 
is the event time and $ is the light intensity per unit solid 

angle (expressed in photons/sr), the amount of light and 
prompt arrival time at the i-th PMT can be predicted as 

fii = ^ei fii Ai (8) 

and 

ff-"^) = + - (9) 

respectively, where is the quantum efficiency of the 
PMT, Qi is the sohd angle subtended by the PMT at a 
distance from the event vertex, Ai is the light trans- 
mission amplitude, and c„ is the effective speed of light 
in the medium. 

The event likelihood is defined as 

£(X) = n /,(0;ai.) n /g(ft;Mi)/t(ii;if'^'\Mz) (10) 

9i=0 9i>0 

where the first product goes over the PMTs that have 
not been hit, while the second product goes over the re- 
maining PMTs that have been hit (i.e., have a non-zero 
recorded charge qi at the registered time U). fq{qi;Hi) 
is the probability to measure a charge qi given an ex- 
pected charge /Uj, and ft{tui''i'^''''^\ P^'i) is the probability 
to measure a time U given a prompt arrival time tY'^'^ and 
predicted charge ^i. These are obtained from MC sim- 
ulations and verified against the physics and calibration 
data. The task of the event reconstruction is to find the 
best possible set of event parameters Xmi„ which max- 
imizes the event likelihood >C(X), or equivalently, mini- 



mizes the negative log-likelihood function 

F(X) = -ln£(X) = -^ln/,(gi;X) 

i 

-Y,\nMU;X)=Fg{X) + Ft{X). (11) 

9t>0 

Note that the event reconstruction can be performed us- 
ing either one or both of the two terms in the expression 
above, Fg for a charge-only reconstruction, or Ft for a 
time-only reconstruction; utilizing both components en- 
hances the accuracy and stability of the algorithm. 

The performance of the Double Chooz reconstruction 
has been evaluated in situ using radioactive sources de- 
ployed at known positions along the 2;-axis in the target 
volume, and off-axis in the guide tubes. The sources 
arc reconstructed with a spatial resolution of 32 cm for 
i37Cs, 24 cm for ^°Co, and 22 cm for ^^Ge. 



C. Muon tagging and reconstruction 

Cosmic muons passing through the detector or the 
nearby rock induce backgrounds which are discussed in 
the next section. A through-going (stopping) muon typi- 
cally deposits 160 MeV (80 MeV) in the IV which triggers 
above about 10 MeV. The IV trigger rate is 46 s~^. All 
muons in the ID are tagged by the IV except some stop- 
ping muons which enter the chimney. Muons which stop 
in the ID and their resulting Michel e can be identified 
by demanding a large energy deposition (roughly a few 
tens of MeV) in the ID. An event is tagged as a muon if 
there is > 5 MeV in the IV or > 30 MeV in the ID. 

Several tracking algorithms have been developed to re- 
construct these muons. IV reconstruction is based on a 
maximum likelihood algorithm utilizing the arrival times 
of the earliest photons to hit each PMT, while ID recon- 
struction utilizes the spatial pattern of hit times. The 
forward wavefront of scintillation light from a relativistic 
track propagates at the Cerenkov angle, thus allowing the 
same algorithm to be used for tracks in the NT, GC, and 
non-scintillating buffer. Using MC and the OV as ref- 
erence, the lateral resolution at the detector center has 
been determined to be 35 cm for ID and 60 cm for IV 
muons. 



D. Light Noise Rejection 

The background known as light noise is caused by a 
sporadic spontaneous flashes of some PMT bases. The 
characteristic signature is light mainly localized to one 
PMT base and spread out in time among the other PMTs 
after many reflections from the detector surfaces. This 
background can be discriminated from physics events 
based on the fact that the detected light is spread less 
homogeneously across the detector for light noise events. 
Light noise is rejected by demanding both a small value of 
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Qmax/Qtot, where Qmax is the maximum charge recorded 
by a single PMT and Qtot is the total ID charge collected 
in a trigger, and large values of rms(istait), which is the 
standard deviation of the distribution of the start time 
(^start) of the first pulse on each PMT. 

E. Energy Reconstruction 

The visible energy (E^ig) provides the absolute calori- 
metric estimation of the energy deposited per trigger. 
Ems is a function of the calibrated PE (total number of 
photoelectrons): 

E,,s = PE"^ip, z, t) X /-(p, z) X /r (t) X (12) 

where PE — pci — qi/gaini{qi). Coordinates in 
the detector are p and z, t is time, m refers to data or 
Monte Carlo (MC) and i refers to each good channel. 
The correction factors /,;, fs and JmcV correspond, re- 
spectively, to the spatial uniformity, time stability and 
EE/MeV calibrations. Four stages of calibration are car- 
ried out to render E^s linear, independent of time and 
position, and consistent between data and MC. Both the 
MC and data are subjected to the same stages of calibra- 
tion. 

The sum over all good channels of the reconstructed 
raw charge (gi, see Section llV A[) from the digitized wave- 
forms is the basis of the energy estimation. Good chan- 
nels are those identified and tagged as well behaved by 
fast online analysis based on waveform information. Only 
a very few channels are sporadically not good and are, 
thus, excluded from the calorimetric estimation. The lim- 
ited sampling of the waveform baseline estimation can 
be biased [2^ leading to a non-linearity at about 1 PE 
charge equivalent. Figure 2] shows the effect for a rep- 
resentative channel. A similar curve is used to calibrate 
the MC. The PE calibrated charge {pei) is defined as 
pei = qi/gaini{qi). One gaini{qi) curve is generated 
upon each power-cycle episode. Due to the average light 
level (^ 230 PE /MeV), the non- linear bias of the single- 
PE calibration can have up to a 10% effect for energies 
below 3 MeV, if not corrected. 

The PE response is position dependent for both MC 
and data. Calibration maps were created such that 
any PE response for any event located at any posi- 
tion {p,z) can be converted into its response as if mea- 
sured at the center of the detector {p = 0,z — 0): 
PE™ = PE"'{p, z) X /™(p, z). The calibration map's cor- 
rection for each point is labeled /™(p, z). Independent 
uniformity calibration maps /™(p, z) are created for data 
and MC, such that the uniformity calibration serves to 
minimize any possible difference in position dependence 
of the data with respect to MC. The capture peak on 
H (2.223 MeV) of neutrons from spallation and antineu- 
trino interactions provides a precise and copious calibra- 
tion source to characterize the response non-uniformity 
over the full volume (both NT and GC). The calibra- 
tion map for data is shown in Figure [5] A similar map 
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FIG. 4. Demonstration of the linear PE calibration for one 
channel. The gain versus charge is shown. The dashed line 
highlights the constant component (linear behavior) of the 
gain observed at large charges. The calibration parametrizes 
this curve to correct the non-linear component (deviation 
from constant) of the gain, making the PE corrected energy 
scale linear to within 2%. 
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FIG. 5. Detector cahbration map, in cylindrical coordinates 
{p,z), as sampled with spallation neutrons capturing in H 
across the ID. Response variations are quantified as the frac- 
tional response with respect to the detector center. Largest 
deviation in NT are up to 5%. A similar map is constructed 
with MC for calibration of its slightly difi^erent response uni- 
formity pattern. 



was measured and applied to MC. A 2D-interpolation 
method was developed to provide a smooth application 
of the calibration map at any point (p, z). 

The detector response stability was found to vary in 
time due to two effects, which are accounted for and cor- 
rected by the term f™{t). First, the detector response 
can change due to variations in readout gain or scin- 
tillator response. This effect has been measured as a 
+2.2% monotonic increase over I year using the response 
of the spallation neutrons capturing on Gd within the 
NT, shown in Figure [H Second, a few readout chan- 
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FIG. 6. Stability of the reconstructed energy as sampled by 
the evolution in response of the spallation neutron H-capture 
after stability calibration. The observed steps correspond to 
power-cycle periods. The systematic uncertainty on the en- 
ergy stability is estimated at 0.61%. 



the non-uniformity systematic uncertainty, and is 0.43%. 
This result is consistent with the analysis of all cali- 
bration sources along the z-axis (NT) and GT (GC). 
The relative instability systematic error, discussed above, 
is 0.61%. Responses are equalized at 2.223 MeV, but 
small data/MC discrepancies in the absolute energy scale 
can still arise from the relative non-linearity across the 
prompt energy spectrum. This possibility was explored 
by using all calibration sources in the energy range 0.7 - 
8 MeV with deployments along the z-axis and GT. Some 
relative non-linearity was observed (< 0.2%/MeV) but 
the pattern diminished when integrated over the full vol- 
ume. A 0.85% variation consistent with this non-linearity 
was measured with the z-axis calibration system, and this 
is used as the systematic error for relative non-linearity 
in Table Hill Consistent results were obtained when sam- 
pling with the same sources along the GT. 



NEUTRINO DATA ANALYSIS 



TABLE III. Energy scale systematic errors. 





Error (%) 


Relative Non-Uniformity 


0.43 


Relative Instability 


0.61 


Relative Non-Linearity 


0.85 


Total 


1.13 



nels varying over time are excluded from the calorimetry 
sum, and the average overall response decreases by 0.3% 
per channel excluded. The MC is stable, so this correc- 
tion is applied only to data. The stability calibration is 
relative to a specific reference time to- Therefore, any 
response PEQ{t) is converted to the equivalent response 
at to, as PEg't^ = PE^it) x The was defined 

as the day of the first Cf source deployment, during Au- 
gust 2011. The remaining instability after calibration is 
shown in Figure |6l as sampled with H-capture from spal- 
lation neutrons, and is used for the stability systematic 
uncertainty estimation. 

The number PEqiq psr MeV is determined by an ab- 
solute energy calibration independently, for the data and 
MC. The response in PEqio for H-capture as deployed 
in the center of the NT is used for the absolute en- 
ergy scale. The absolute energy scales are found to be 
229.9 PEQtjMeY and 227.7 PE^tJUeY, respectively, 
for the data and MC, demonstrating agreement within 
1% prior to this calibration stage. 

Discrepancies in response between the MC and data, 
after calibration, are used to estimate these uncertain- 
ties within the prompt energy range and the NT vol- 
ume. Table IIIII summarizes the systematic uncertainty 
in terms of the remaining non-uniformity, instability and 
non-linearity. The relative non-uniformity systematic un- 
certainty was estimated from the calibration maps using 
neutrons capturing on Gd, after full calibration. The 
rms deviation of the relative difference between the data 
and MC calibration maps is used as the estimator of 



A. Ve Candidate Selection 

The candidate selection procedure starts in a similar 
way as [7(1 • Events with an energy below 0.5 MeV, where 
the trigger efficiency is not 100%, or identified as light 
noise (Qmax/Qtot > 0.09 or rms(istait) > 40 ns) are dis- 
carded. Triggers within a 1 ms window following a tagged 
muon are also rejected (see Section HV C|) . in order to re- 
duce the correlated and cosmogenic backgrounds. The 
effective veto time is 4.4% of the total run time. Defin- 
ing AT = ^delayed " ^prompt , further selection consists of 
4 cuts: 

1. time difference between consecutive triggers 
(prompt and delayed): 2 /is < AT < 100 /iS, as 
shown in Figure [3 where the lower cut reduces 
correlated backgrounds and the upper cut is 
determined by the approximately 30 /xs capture 
time on Gd; 

2. prompt trigger: 0.7 MeV < £;prompt < 12.2 MeV, 
as illustrated in Figure [5] 

3. delayed trigger: 6.0 MeV < ^^dciaycd < 12.0 MeV 
(Figured]) and Qmax/Qtot <0.055; 

4. multiplicity: no additional triggers from 100 /iS pre- 
ceding the prompt signal to 400 /iS after it, with the 
goal of reducing the correlated background. 

The IBD efhciencies for these cuts are listed in Table HVl 

A preliminary sample of 9021 candidates is obtained by 
applying selections (1-4). In order to reduce the back- 
ground contamination in the sample, candidates are re- 
jected according to two extra cuts not used in Q- First, 
candidates within a 0.5 s window after a high energy 
muon crossing the ID {Efj_ > 600 MeV) are tagged as 
cosmogenic isotope events and rejected, increasing the 
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Cut 


Efficiency % 


-f/prompt 
-^delayed 

A T 
Multiplicity 
Muon veto 
Outer Veto 


100.0 ± 0.0 

94.1 ± 0.6 

96.2 ± 0.5 
99.5 ± 0.0 

90.8 ± 0.0 

99.9 ± 0.0 



TABLE IV. Cuts used in the event selection and their ef- 
ficiency for IBD events. The OV was working for the last 
68.9% of the data. 



I I anti-v MC 
• Data 




AT (us) 

FIG. 7. Time difference between prompt and de- 
layed triggers. Black dots and solid histogram 
show data and MC results, respectively. 



effective veto time to 9.2%. Second, candidates whose 
prompt signal is coincident with an OV trigger are also 
excluded as correlated background. Applying the above 
vetoes yields 8249 candidates or a rate of 36.2 ± 0.4 
events/day, uniformly distributed within the target, for 
an analysis livetime of 227.93 days. This rate is lower 
than the one presented in Q due to a longer data taking 
period with one reactor being off, as well as to the new 
cuts reducing the background contamination. Following 
the same selection procedure on the MC sample yields 
8439.6 expected events in the absence of oscillation. 



B. Accidental Background 

The main source of accidental coincidences is the ran- 
dom association of a prompt trigger from natural radioac- 
tivity and a later neutron-like candidate. This back- 
ground is estimated by applying the neutrino selection 
cuts described in Section IVAI but using coincidence win- 
dows shifted by 1 s in order to remove correlations in the 
time scale of n-captures in H and Gd. The statistics of 




prompt E (MeV) 

FIG. 8. Delayed energy versus prompt energy for 
time-correlated triggers. Vertical and horizontal 
dashed lines show the cuts applied for i/^ candi- 
dates selection. 



the sample is enhanced by using 198 windows each shifted 
from the previous one by 500 /iS. The radioactivity rate 
between 0.7 and 12.2 MeV is 8.2 s~^, while the singles 
rate in 6-12 MeV energy region is 18 h~^. Both rates 
are quite stable along the data taking period. Finally, the 
accidental background rate is found to be 0.261 ± 0.002 
events per day. The reproducibility of our result and any 
possible systematic effect are studied by repeating the 
procedure 30 times, i.e., taking 30 times 198 consecutive 
time windows. The dispersion of these 30 measurements 
is consistent with only statistical error, so, no systematic 
deviation is found. 

Figure |9] shows the accidental prompt spectrum and 
the energy distribution for natural radioactivity scaled 
to the number of accidental events; the agreement is 
excellent. The distribution is peaked at low energies 
below 3 MeV. The remaining light noise is included in 
the accidental background sample. Using the correlation 
between both variables Qmax/Qtot and rms(tstart), its 
contribution to the accidental sample is estimated to be 
lower than 1%. 



C. Cosmogenic Isotopes Background 

The radioisotopes ^He and ^Li are products of spalla- 
tion processes on ^^C induced by cosmic muons crossing 
the scintillator volume. The /3n-decays of these isotopes 
constitute a background for the antineutrino search. I3n- 
emitters can be identified from the time- and space- cor- 
relation to their parent muon. Due to their relatively 
long lifetimes (^Li: r ~ 257 ms, ®He: t = 172 ms), an 
event-by-event discrimination is not possible. For the 
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FIG. 9. The accidental prompt spectrum (black circles) su- 
perimposed to the radioactivity energy distribution measured 
in Double Chooz scaled to the same number of entries (red 
line). 




10 12 
Prompt Energy [MeV] 

FIG. 10. The prompt j3 spectrum of the /3n-emitters ^He and 
^Li from data (black squares) and Monte Carlo (red line), 
assuming 



muon rates in our detector, vetoing for several isotope 
lifetimes after each muon would lead to an unacceptably 
large loss in exposure. Instead, the rate is determined 
by an exponential fit to the l^t^^ = — profile of all 
possible muon-IBD candidate-pairs. 

The analysis is performed for three visible energy E"^'" 
ranges that characterize subsamples of parent muons by 
their energy deposition, not corrected for energy non- 
linearities, in the ID: 

1. Showering muons crossing the target value are se- 
lected by E^^'^ > 600 MeV and feature an increased 
probability to produce cosmogenic isotopes. The 
Aip^-fit returns a precise result of 0.95 ±0.11 
events/day for the /3n-emitter rate. 

2. In the range from 275 to 600 MeV, muons 
crossing GC and target still give a sizable con- 
tribution to isotope production of 1.08 ± 0.44 
events/day. To obtain this result from a At^i, fit, 
the sample of muon-IBD pairs has to be cleaned 
by a spatial cut on the distance of closest ap- 
proach from the muon to the IBD candidate of 
df^i, < 80 cm to remove the majority of uncorre- 
lated pairs. The corresponding cut efficiency is de- 
termined from the lateral distance profile obtained 
for E^^'" > 600 MeV. The approach is validated by 
a comparative study of cosmic neutrons that show 
an almost congruent profile with very little depen- 
dence on E^;' above 275 MeV. 

3. The Cut E^"^" < 275 MeV selects muons crossing 
only the buffer volume or the rim of the GC. For 
this sample, no production of /3n-emitters inside 
the target volume is observed. An upper limit of < 
0.3 events/day can be established based on a Atf^^, 
fit for dfj_i, < 80 cm. Again, the lateral distribution 
of cosmic neutrons has been used for determining 
the cut efficiency. 



The overall rate of /3n-decays found is 2.05_g'52 
events/day. The result of a similar analysis based on 
the IV muon tracking agrees within the uncertainty. 

Accidental coincidences containing the /3-decay of the 
isotope ^^B either as prompt or as delayed event feature 
a time correlation to the parent muons producing the 
^^B. In the At^^y profile, these events are represented by 
a decay function with t{^'^B) — 29 ms. However, these 
events were removed very efficiently from the data set 
used for ^Li analysis by imposing a maximum distance 
cut of 90 cm between prompt and delayed events, intro- 
ducing a negligible inefficiency of ^1%. 

The correlation of a cosmogenic isotope to the show- 
ering muons has been exploited to impose a partial 
veto of this background for the final fit analysis. Ve- 
toing all IBD candidates within 0.5 s following a muon 
of El^' > 600 MeV, 0.89 ± 0.10 events d~i of ;3n-decays 
are removed from the data sample. The residual cosmo- 
genic isotope background rate has been determined to 
1.25 ± 0.54 events/day. 

Finally, the correlation of parent muons and l3n- 
emitters has been used to extract the prompt /? spectrum 
from the data. Figure[TO] shows a sample spectrum ob- 
tained for Efj_ > 620 MeV, a distance cut of 0.7 m and 
a At^^ cut of 600 ms. The contamination of the sample 
by random coincidences has been statistically subtracted. 
Good agreement is found for the MC spectrum used in 
the final fit analysis. 

D. Fast Neutrons and Stopping Muon Background 

Most correlated backgrounds are rejected by the 1 ms 
veto time after each tagged muon. The remaining events 
arise from cosmogenic events whose parent muon either 
misses the detector or deposits an energy low enough to 
escape the muon tagging. Two contributions have been 
found: fast neutrons (FN) and stopping muons (SM). 
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FN are created by muons in the inactive regions sur- 
rounding the detector. Their large interaction length al- 
lows them to cross the detector and capture in the ID, 
causing both a prompt trigger by recoil protons and a 
delayed trigger by capture on Gd. An approximately flat 
prompt energy spectrum is expected; a slope could be 
introduced by acceptance and scintillator quenching ef- 
fects. The time and spatial correlation distribution of FN 
are indistinguishable from those of events. 

The selected SM arise from muons entering through 
the chimney, stopping in the top of the ID, and eventu- 
ally decaying. The short muon track mimics the prompt 
event, and the decay Michel electron mimics the delayed 
event. SM candidates are localized in space in the top 
of the ID under the chimney, and have a prompt-delayed 
time distribution following the 2.2 /is muon lifetime. 

The correlated background has been studied by ex- 
tending the selection on E'prompt up to 30 MeV. No IBD 
events are expected in the interval 12 MeV < ii'prompt < 
30 MeV. FN and SM candidates were separated via 
their different correlation time distributions. A 97l;g% 
pure sample of FN is obtained for AT > 10 /is, and a 
(88 ± 7)% pure sample of SM is obtained for AT < 10 ^is. 
These samples of FN and SM can be used to estimate 
their rate. The observed prompt energy spectrum is con- 
sistent with a flat continuum between 12 and 30 MeV, 
which extrapolated to the IBD selection window provides 
a first estimation of the correlated background rate of 
« 0.75 events/day. The accuracy of this estimate de- 
pends on the validity of the extrapolation of the spectral 
shape. Next we describe a measurement of the FN and 
SM spectral shapes including the IBD region, obtained 
by using the IV and OV to tag samples of FN and SM. 

The DAQ reads out the IV upon any ID trigger, lower- 
ing the IV detection threshold to ^ 1 MeV, and making 
the IV sensitive to FN via the detection of proton recoils 
and captures on H. The IV-tagging is implemented by 
demanding at least 2 IV PMT hits leading to (33 ± 5)% 
tagging efficiency with no contribution by single PMT 
energy depositions. There is a very low probability of 
accidental IV tagging due to any IV energy deposition in 
the 256 ns coincident readout window. 

The OV-tagging, when available, is especially sensitive 
to SMs since the muon is often detected. (41 ± 23)% of 
the FN and SM candidates in the 12-to-30-MeV window 
are tagged by the OV, of which (74± 12)% are SM. OV- 
tagging has an accidental rate = 0.06% of the neutrino 
sample and can be used to veto events caused by muons. 

Several FN and SM analyses were performed using dif- 
ferent combinations of IV and OV tagging. The main 
analysis for the FN estimation relies on IV-tagging of the 
prompt triggers with OV veto applied for the IBD selec- 
tion. Two sources of backgrounds on the tagged FN sam- 
ple were identified and rejected. The first source is the 
combination of natural radioactivity in the IV in an acci- 
dental coincidence with a genuine IBD, and was reduced 
to 12% by imposing a time coincidence between the ID 
and IV energy depositions. The second source, a Comp- 




Energy (MeV) 

FIG. 11. FN and SM combined spectral model best fit (solid 
red) with ± la (dashed red) , energy distribution of tagged FN 
and SM population (gray histogram) and IBD spectrum. 



ton scattering in both the IV and ID in an accidental 
coincidence with a Gd-capture, was reduced to 2% by im- 
posing a cut on the spatial distance between the prompt 
and delayed candidate in the ID. The purity of the IV- 
tagged FN sample was 86%. The remaining background 
was measured in an off-time window and subtracted, thus 
minimizing distortions to the energy spectrum. The FN 
spectral shape was found to be in agreement with a linear 
model with a small positive slope. The measured total 
FN rate was (0.30 ± 0.14) events/day, including system- 
atic uncertainties from the AT-based FN-SM separation, 
the IV-tagging efficiency, and background subtractions. 

Since there is no correlation between the SM prompt 
energy and the delayed energy deposit of the Michel elec- 
tron, a pure sample of SM was obtained by selecting 
20 MeV < T^dciayod < 60 MeV. The spectral shape of 
SM prompt energy was found to be in agreement with 
a linear model with a small negative slope. The total 
SM rate was measured to be (0.34 ± 0.18) events/day, 
including systematic uncertainties. 

Since the spectral shapes for both FN and SM are lin- 
ear, a combined analysis was performed to obtain the 
total spectrum shown in Fig. [TT] and the total rate esti- 
mation (0.67±0.20) events/day summarized in Table fVl 
Consistent results were obtained from different analysis 
techniques, which included IV- and OV-tagging without 
OV- vetoing. The OV veto reduces the rate of correlated 
backgrounds by about 30%. 

E. Background Measurements 

There are four ways that can be utilized to estimate 
backgrounds. Each independent background component 
can be measured by isolating samples and subtracting 
possible correlations. This is described for each compo- 
nent in Sections II V D 1 IVC] and IVDl Second, we can mea- 
sure each independent background component including 
spectral information when fitting for ^13 oscillations as is 
done in Section IVII Third, the total background rate is 
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measured by comparing the observed and expected rates 
as a function of reactor power. Fourth, we can use the 
both-reactor-off data to measure both the rate and spec- 
trum. 

The latter two methods are used currently as cross- 
checks for the background measurements due to low 
statistics and are described here. The measured daily 
rate of IBD candidates as a function of the no-oscillation 
expected rate for different reactor power conditions is 
shown in Figure [121 The extrapolation to zero reactor 
power of the fit to the data yields 2.9 ± 1.1 events per 
day, in excellent agreement with our background esti- 
mate. The overall rate of correlated background events 
that pass the IBD cuts is independently verified by ana- 
lyzing 22.5 hours of both-reactors-off data. The expected 
neutrino signal is < 0.3 residual Ve events. Three events 
passed the first 4 cuts in Section fV Al Two events with 
prompt energies of 4.8 MeV and 9.4 MeV were associated 
within 30 cm and 240 ms with the closest energetic muon, 
and are thus likely to be associated with ^Li. Indeed, the 
second candidate is rejected by the showering muon veto. 
The third candidate at a prompt energy of 0.8 MeV fea- 
tures 3.5 m distance between prompt and delayed events 
and is therefore most likely a random coincidence. Im- 
mediately following the data set used in this paper, we 
obtained a larger data set with both-reactors-off. That 
will be the subject of a separate paper [59j . 
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FIG. 12. Daily number of />e candidates as a func- 
tion of the expected number of Pg. The dashed 
line shows the fit to the data, along with the 90% 
C.L. band. The dotted line shows the expectation 
in the no-oscillation scenario. 



F. Neutron Detection Efficiency 

Calibration data taken with the ^^^Cf source were used 
to check the Monte Carlo prediction for any biases in 
the neutron selection criteria and estimate their contri- 
butions to the systematic uncertainty. 

The fraction of neutron captures on gadolinium is eval- 
uated to be 86.5% near the center of the target, 1.5% 
lower than the fraction predicted by simulation. There- 
fore the Monte Carlo simulation for the prediction of the 
number of Pe events is reduced by factor of 0.985. Af- 
ter the prediction of the fraction of neutron captures on 
gadolinium is scaled to the data, the prediction repro- 
duces the data to within 0.3% under variation of selection 
criteria. 

The ^^^Cf is also used to check the neutron capture 
time, AT. The time difference between the prompt event 
and neutron capture signal for the californium calibration 
data is shown in Figure 1131 The simulation reproduces 
the efficiency (96.2%) of the Aig+„ cut with an uncer- 
tainty of 0.5% augmented with sources deployed through 
the NT and GC. 

The efficiency for Gd capture events with visible energy 
greater than 4 MeV to pass the 6 MeV cut is estimated 
to be 94.1%. Averaged over the NT, the fraction of neu- 
tron captures on Gd accepted by the 6.0 MeV cut is in 
agreement with calibration data to within 0.7%. 

The Monte Carlo simulation indicates that the number 
of IBD events occurring in the GC with the neutron cap- 
tured in the NT (spill-in) slightly exceeds the number of 
events occurring in the target with the neutron escaping 
to the gamma catcher (spill-out), by 1.35% ± 0.04%(stat) 
± 0.30%(sys). The spill-in/out effect is already included 
in the simulation and therefore no correction for this is 
needed. The uncertainty of 0.3% assigned to the net spill- 
in/out current was quantified by varying the parameters 
affecting the process, such as gadolinium concentration 
in the target scintillator and hydrogen fraction in the 
gamma-catcher fluid within its tolerances. Moreover the 
parameter variation was performed with multiple Monte 
Carlo models at low neutron energies. 

VI. OSCILLATION ANALYSIS 

The oscillation analysis is based on a combined fit to 
antineutrino rate and spectral shape. IBD candidates are 
selected as described in Section FV Al The data are com- 
pared to the Monte Carlo signal and background events 
from high-statistics samples. The same selections are ap- 
plied to both signal and background, with corrections 
made to Monte Carlo only when necessary to match de- 
tector performance metrics. 

The oscillation analysis begins by separating the data 
into 18 variably-sized bins between 0.7 and 12.2 MeV. 
Two integration periods are used in the fit to help sepa- 
rate background and signal flux. One set contains data 
periods where one reactor is operating at less than 20% of 
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AT(ns) 

FIG. 13. Time difference between prompt and delayed events 
with ^^^Cf at the detector center. The prompt time is deter- 
mined by 7-30 MeV gamma ray. 

its nominal thermal power, according to power data pro- 
vided by EDF, while the other set contains data from all 
other times, typically when both reactors are running. 
All data end up in one of the two integration periods. 
Here, we denote the number of observed IBD candidates 
in each of the bins as iV^ , where i runs over the combined 
36 bins of both integration periods. The use of multiple 
periods of data integration takes advantage of the differ- 
ent signal/background ratios in each period, as the signal 
rate varies with reactor power while the backgrounds re- 
main constant in time. This technique adds information 
about background behavior to the fit. The distribution 
of IBD candidates between the two integration periods is 
given in Table |Vl 

A prediction of the observed number of signal and 
background events is constructed for each energy bin, fol- 
lowing the same integration period division as the data: 

Reactors Bkgnds. 

Nr'= E E (13) 

R=l,2 b 

where N^'^ — Pj-pg Dg)N^^^'^; Pfj^^fj^ is the neutrino 
survival probability from the well-known oscillation for- 
mula and ]\[°^^'^ is given by Equation 21 The index b 
runs over the three backgrounds: cosmogenic isotope; 
correlated; and accidental. The index R runs over the 
two reactors, Chooz Bl and B2. 

Background populations were calculated based on the 
measured rates and the livetime of the detector during 
each integration period. Details on the signal prediction 
normalization can be found in Sec. IIIIDI Predicted pop- 
ulations for both null-oscillation signal and backgrounds 
may be found in Table El 

Systematic and statistical uncertainties are propagated 
to the fit by the use of a covariance matrix Mij in order 
to properly account for correlations between energy bins. 





Reactors 


One Reactor 


Total 




Both On 


Pth < 20% 




Livetime [days] 


139.27 


88.66 


227.93 


IBD Candidates 


6088 


2161 


8249 


!^ Reactor Bl 


2910.9 


774.6 


3685.5 


u Reactor B2 


3422.4 


1331.7 


4754.1 


Cosmogenic Isotope 


174.1 


110.8 


284.9 


Correlated FN & SM 


93.3 


59.4 


152.7 


Accidentals 


36.4 


23.1 


59.5 


Total Prediction 


6637.1 


2299.7 


8936.8 



TABLE V. Summary of observed IBD candidates, with cor- 
responding signal and background predictions for each inte- 
gration period, before any oscillation fit results have been 
applied. 



The sources of uncertainty A are listed in Table I VI I 



Bkgnds. 

= M^^^- + Mf^'- + Af^*"*- + M^f- + ^ M,'^. (14) 

b 



Each term = cov(iVP"'^, iVP"'^)^ on the right-hand 
side of Equation [T31 represents the covariance of Nf'^'"^ 
and NY°'^ due to uncertainty A. The normalization un- 
certainty associated with each of the matrix contribu- 
tions may be found from the sum of each matrix: these 
are summarized in Table IVTl Many sources of uncertainty 
contain spectral shape components which do not directly 
contribute to the normalization error, but do provide for 
correlated uncertainties between the energy bins. The 
signal covariance matrix M^^^' is calculated taking into 
account knowledge about the predicted neutrino spectra. 
The ^Li matrix contribution contains spectral shape un- 
certainties estimated using different Monte Carlo event 
generation parameters, as described in Sec. IIII El The 
slope of the FN/SM spectrum is allowed to vary from 
a nearly-flat spectrum following the measurements de- 
scribed in Section IV Dl Since accidental background un- 
certainties are measured to a high precision from many 
off-time windows, they are included as a diagonal covari- 
ance matrix. 

The elements of the covariance matrix contributions 
are recalculated as a function of the oscillation and other 
parameters (see below) at each step of the minimization. 
This maintains the fractional systematic uncertainties as 
the bin populations vary from the changes in the oscilla- 
tion and fit parameters. 

A fit of the binned signal and background data to 
a two-neutrino oscillation hypothesis was performed by 
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Source 


Uncertainty [%] 


Reactor Flux 


1.67% 


Detector Response 


0.32% 


Statistics 


1.06% 


Efficiency 


0.95% 


Cosmogenic Isotope Background 


1.38% 


FN/SM 


0.51% 


Accidental Background 


0.01% 


Total 


2.66% 



TABLE VI. Summary of signal and background normalization 
uncertainties in tliis analysis relative to the total prediction. 



minimizing a standard function: 



36 



i-.j 



{as - if , (Amii - {Aml^] 



MINOSV 



(15) 



^ MINOS 



The use of energy spectrum information in this analy- 
sis allows additional information on background rates to 
be gained from the fit, in particular because of the small 
number of IBD events between 8 and 12 MeV. The two fit 
parameters epN/SM and egLj are allowed to vary as part 
of the fit, and they scale the rates of the two backgrounds 
(correlated and cosmogenic isotope). The rate of acci- 
dentals is not allowed to vary since its initial uncertainty 
is precisely determined by the measurement method de- 
scribed in Sec. IV Bl The energy scale for predicted signal 
and ^Li events is allowed to vary linearly according to the 
a_B parameter with an uncertainty — 1-13%. A final 
parameter constrains the mass splitting Am§]^ using the 
MINOS measurement ^ of Ato|i = (2.32±0.12) x lO'^ 
eV^, where we have symmetrized the error. This error in- 
cludes the uncertainty introduced by relating the effective 
mass-squared difference observed in a disappearance 
experiment to the one relevant for reactor experiments, 
and the ambiguity due to the type of the neutrino mass 
hierarchy, see e.g. 



Uncertainties for these parame- 
ters, (JpN/SMi <^^Li-: and (TminoSi are listed as the initial 
values in Table EUl 

The best-fit gives sin^ 26li3 = 0.109 ± 0.030 (stat.) ± 
0.025 (syst) at Ato|i = 2.32 x 10"^ gy^, with a 
xVNDF = 42.1/35. We used the MINOS measured 
Am§^ value as a constraint for our ^13 measurement, but 
a two parameter fit without the MINOS Am§i in the re- 
gion Amg]^ < 0.01 eV^ gives a Am|]^ value of 2.7 ± 
1.9 X 10"^ eV^ which is fully consistent with MINOS. 
The fit gives sin^ 26*13 = 0.093 ± 0.078 which is consistent 
with our fit for 6*13 using MINOS. 



Table IVIII gives the resulting values of the fit parame- 
ters and their uncertainties. Comparing the values with 
the ones used as input to the fit in Table IVl we conclude 
that the background rate and uncertainties are further 
constrained in the fit, as well as the energy scale. 

The final measured spectrum and the best-fit spectrum 
are shown in Figure[T3]for the new and old data sets, and 
for both together in Figure [121 



Fit Parameter 


Initial Value 


Best-Fit Value 


'■'Li Bkg. e9Li 
FN/SM Bkg. epM/SM 
Energy Scale q_b 
Amii (10-^ eV^) 


(1.25 ± 0.54) d"" 
(0.67 ± 0.20) d"^ 
1.000 ± 0.011 
2.32 ± 0.12 


(1.00 ± 0.29) d"' 
(0.64 ± 0.13) d"^ 
0.986 ± 0.007 
2.32 ± 0.12 



TABLE VII. Parameters in the oscillation fit. Initial val- 
ues are determined by measurements of background rates or 
detector calibration data. Best-fit values are outputs of the 
minimization procedure. 



An analysis comparing only the total observed number 
of IBD candidates in each integration period to the expec- 
tations produces a best-fit of sin^ 2^13 = 0.170 ± 0.052 at 
xVNDF = 0.50/1. The compatibility probability for the 
rate-only and rate-|-shape measurements is about 30% 
depending on how the correlated errors are handled be- 
tween the two measurements. 

A re-processing of the data set used for the first Dou- 
ble Chooz publication [7] was performed using the cur- 
rent analysis techniques. A fit using only a single in- 
tegration period yielded a best-fit value of sin^ 2^13 = 
0.0744 ± 0.046 with xV^DF = 18.3/17. An anal- 
ysis of only the data taken since the first publication 
yielded a best-fit of sin^ 26'i3 = 0.143 ± 0.043 with 
xVNDF = 9.54/17. The data and best-fit spectra for 
each of these cases is shown in Figure [121 

Our predicted fission cross section is 5.723 ± 0.096 
X lO^'*'^ cm^/fission using the Bugey4 anchoring measure- 
ment and corresponding to the values of ak in Table HI 
The background subtracted reactor antineutrino event 
rate is 7751.9 events, corresponding to 91.85% of the rate 
expected in the absence of oscillations. Our measured fis- 
sion cross section is 5.257 ± 0.056 (stat.) ± 0.105 (syst) 
xlO~^'^ cm^ /fission. 

A further cross-check of the analysis was carried out 
by imposing cuts to eliminate the vast majority of the 
cosmogenic isotope background at the cost of reduced 
livetime. The best-fit case of this analysis was found at 
sin^ 26li3 = 0.109 ± 0.044 and Am§i = 2.32 x lO^^ eV^, 
in good agreement with the standard analysis. 

Confidence intervals for the standard analysis were de- 
termined using a frequentist technique [g^l- This ap- 
proach accommodates the fact that the true x^ distribu- 
tions may not be Gaussian and is useful for calculating 
the probability of excluding the no-oscillation hypoth- 
esis. This study compared the data to 10,000 simula- 
tions generated at each of 21 test points in the range 
< sin^ 26*13 < 0.25. A Ax^ statistic, equal to the dif- 
ference between the x^ at the test point and the x^ at 
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FIG. 14. Measured prompt energy spectrum for each integration period (data points) superimposed on the expected prompt 
energy spectrum, including backgrounds (green region), for the no-oscillation (blue dotted curve) and best-fit (red solid curve) 
at sin^ 2613 = 0.109 and Amli — 2.32 X 10 eV^. Inset: stacked spectra of backgrounds. Bottom: differences between data 
and no-oscillation prediction (data points), and differences between best fit prediction and no-oscillation prediction (red curve). 
The orange band represents the systematic uncertainties on the best-fit prediction. 




FIG. 15. Sum of both integration periods plotted in the same 
manner as Figure [HI 



the best fit, was used to determine the region in sin^ 20i3 



where the Ax^ of the data was within the given confi- 
dence probability. The allowed region at 68% (90%) CL 
is 0.067 (0.043) < sin^ 26*13 < 0.15 (0.18). An analogous 
technique shows that the data excludes the no-oscillation 
hypothesis at 99.8% (2.9a). 



VII. CONCLUSION 

A comparison of this analysis result to other recent 
sin^ 2013 measurements by other experiments is shown in 
Figure 1171 The figure shows published results, though 
we note that new results from Daya Bay, MINOS and 
T2K have been shown at conferences but are not yet 
published [63i]. The values for sin^ 26*13 from the various 
experiments are in excellent agreement with the results 
reported here. However this result is unique in its incor- 
poration of energy dependence in the analysis. 

Double Chooz has found evidence for a non-zero value 
of ^13 from the rate and energy spectrum of reactor neu- 
trino candidates at a distance of 1050 m from two re- 
actors. It is the first evidence for this parameter using 
the energy spectrum from reactor neutrinos, rather than 
simply their rate. We find a best fit value and Icr error to 
be sin^ 26*13 = 0.109 ± 0.030 (stat) ± 0.025 (syst). The 
data is inconsistent with the assumption that oscillations 



20 



5- 1000 

« 800 
c 

iS 700 

600 
500 
400 
300 
200 



100 3FI 



1.4 
1.2 
1.0 
0.8 
0.6 
50 

0, 
-50 
-100 



Double Chooz 201 1 Data Set 
No Oscillation, Best-fit Bacitgrounds 
Best Fit: sin^(2e^g) = 0.074 
at Aitlg, = 0.00232 eV^* (x^/d.o.f. = 18.3/17) 
Summed Backgrounds (see insets) 
Lithium-9 

Fast n and Stopping \i 
Accidentals 




"Ttr" — — ^z 

Energy [MeV] 



(a)First Publication Data Set 



J- 1000 
S 000 
^ 800 





m 700 
600 
500 
400 
300 
200 



Double Chooz Data Since Nov. 2011 

No Oscillation, Best-fit Backgrounds 

Best Fit: sin^(2e^g) = 0.143 

at Am^, = 0.00232 eV^ (x^/d.o.f. = 9.5/17) 

Summed Backgrounds (see insets) 

Lithium-g 

Fast n and Stopping p 
Accidentals 





(b)Data Since First Publication 



"Ttr" — — ^Z 
Energy [MeV] 



FIG. 16. Data and best-fit spectrum from applying current 
analysis techniques to the data set used to produce the first 
Double Chooz publication (a), and data taken since that pub- 
lication (b), plotted in the same manner as Figure [141 



are absent with a CL of 99.8% CL (2.9cr). 
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